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ABSTRACT 



Hard X-ray INTEGRAL observations of SScarried out 
during 2003-2005 years are presented. Analysis of pre- 
cessional and orbital variability is presented. The width 
of X-ray eclipse in the 25 — 50 keV range at the preces- 
sional phase tp = 0.1 (accretion disk is open to observer) 
is higher than that in the Ginga 18.4 — 27.6 keV range. 
This fact suggests existance the presence of hot extended 
corona around the supercritical accretion disk. Spectrum 
of hard X-rays in the range 10 — 200 keV does not change 
with the precessional phase which also suggests that hard 
X-ray flux is generated in the hot extended corona around 
the accretion disk. The parameters of this hot corona are: 
kT = 23 - 25 keV, r = 1.8 - 2.8. Mass ratio estimated 
from the analysis of the ingress part of the eclipse light 
curve is in the range q — m x /m v = 0.3 — 0.5. 



1. INTRODUCTION 



SSis a massive eclipsing X-ray binary with precessing su- 
percritical accretion disk and relativistic jets. Narrow col- 
limated relativistic jets (v = 0.26c) are precessing with 
the period P prec = 162 d .5. Orbital period is 13 d .08 (H). 
There are problems with optical classification of the de- 
termination of its radial velocity curve ( |2|. @], @|, @|, 
((J). This unique X-ray binary at an advanced evolution- 
ary stage has been investigated in optical, radio and X-ray 
ranges (see review of Fabrika |7] and references therein). 

First INTEGRAL observations of SSin hard X-rays gave a 
surprise: SSis a hard X-ray source with emission clearly 
detected up to 100 keV. We concluded that SSis a super- 
critical microquasar with hard X-ray spectrum (| 8;]). SS- 
was observed by INTEGRAL in A01-A03 from 2003 to 
2005. 



2. OBSERVATIONS AND DATA REDUCTION 



In this work all publically available data on SSobtained 
by INTEGRAL from March, 2003 to November, 2004 
(AO- 1,2) and the results of our pointed observations at 
October 12-22, 2005 (AO-3) are presented. 

The data reduction was performed by the Offline 
Science Analysis (OS A) software version 5.1, devel- 
oped by the INTEGRAL science data center (ISDC, 
http://isdc.unige.ch QJ). 

The light curves where reconstructed by processing all 
the available observations, obtained by the ISGRI/IBIS 
detector in the 25 — 50 keV energy band. The light curve 
points where taken from the mosaic images integrated by 
each 10 ScW (Science Windows) with a standard back- 
ground approximation procedure. Thus, a typical exposi- 
tion time for each light curve point is ~ 20 ksec. 

The variability analysis is based on the precessional and 
orbital ephemeris of SSby |7]. The orbital minima: 

JD Mini He = 2450023.62+ 13.08211 * E 

and the maximal emission lines separation T3 moments: 

JD T3 = 2443507.47+ 162.375 * El, 

where E and El are the orbital and precession cycle 
numbers respectively. 

During the 3 years of INTEGRAL observations (2003- 
2005) cover a wide range of the precession phases of SS. 
The three observation sets of AO-1 are performed at ip = 
0.72-0.88 (incontiunosly), 0.95-0.1 (continuosly); AO- 
2: ij) = 0.0 -0.15 (incontinuously), 0.22 -0.23 (continu- 
ously), 0.33-0.35 (continuously); AO-3: ip = 0.5-0.56. 
The time distribution of these observations is displayed in 
FigH The T3 moments are indicated by squares and the 
closest orbital minima - by triangles. 
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Figure 1. Observations of SSin 2003 — 2005. The squares 
indicate the precessional face-on (T3) moments, and the 
triangles indicate the orbital minima. The units of the 
horizontal axis are IJD = JD — 2451544.5. 



3. PRECESSIONAL VARIABILITY 



The 25 — 50 keV light curve composed from all the avail- 
able ISGRI/IBIS observations of SSfolded with the pre- 
cessional period P prec = 162 d .5 is displayed in Fig.|2] 

The face-on disk position (T3) defined as the maximim 
of the optical emission lines separation is taken as the 
moment of zero precession phase. The upper and lower 
panels of Fig. [2] display the X-ray flux (mCrabs) out of 
eclipse of the X-ray source by the normal star and in the 
middle of the eclipse, respectively. 

Fig. |3 shows that the flux out of eclipse is varies from 
3 mCrab (crossover phase) to 18 mCrab (T3 phase). The 
flux in the middle of eclipse (primary minimum) has no 
significant variablity, staying at average near 2.9 mCrab. 



4. ORBITAL ECLIPSE VARIABILITY 



Fig. [3] displays the shape of the eclipses of SSX-ray 
source by the normal component in the energy band 
25 — 50 keV observed in different phases ip of the pe- 
riod of precession. The eclipse at ip = 0.1 is obtained 
by folding of all the available light curves in the phases 
4> = 0.95 - 0.15 (JD 2452770-2453270). 

The eclipse depth is maximal at ip = 0.1 and decreases 
with the phase of precession. This seen clearly for in- 
dividual eclipse light curves for ip — 0.22, 0.34 (first 
crossover), ip — 0.5, 0.72 (some later than the second 
crossover), and 0.88 (see Fig.|3jl. The dotted line at the 
figure indicate the moments of the primary minimum. 
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Figure 2. Precessional light curves ofSS 



The good correspondence between X-ray minima with 
the optical ephemeris for the precession phases close to 
T3 as well as the delay of the X-ray light curve after the 
optical one at ip ~ 0.5 should be specially noted. 
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Figure 3. Orbital light curves of SSin different phases ip 
of precession. 



5. ECLIPSE IN THE PHASE ip = 0. 1 

Fig. ^displays the eclipse light curves in the band 25 — 
50 keV obtained by INTEGRAL in comparison with the 
18.4 - 27.6 keV light curves obtained by Ginga (floh. 
denoted below as hard and soft bands. 

Evidently, the eclipse width in the hard band is greater 
than in the soft one. This seems unexpected because 



Table 1. Parameters of the SS433 spectrum for CompPS model 



Model 


kT (keV) 


T v 


H/R 


Normalization 


Reduced x 2 


hemisphere 


24.9 ±7.7 


1.8 ±0.6 




13.3 ±3.8 


1.08 (5 d.o.f.) 


slab 


23.2 ±6 


2.8 ±0.6 


0.5 


12.7 ±4.5 


1.06 (5 d.o.f.) 
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GINGA LIGHT CURVE SS433 V = 0,1 [18.4-27.6 keV] 
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INTEGRAL CONSOLIDATE LIGHT CURVE SS433 14 = 0,1 [25-50 keV] 
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Figure 5. The IBIS/ISGRI 10 - 200 keV high and low 
state SSspectra (upper and lower points, respectively). 



Figure 4. The eclipse shape at ip = 0.1 



such behaviour is inverse to that in regular X-ray binaries 
where the X-ray eclipse width decreases in harder band. 

The ascending part of the eclipse in the Ginga data is dis- 
turbed, probably, by the light absorption by the accretion 
flows. Some distortions at the same phases can be noticed 
also in the hard band data. 

The phase of the excess point inside the eclipse is coinci- 
dent with the local maximum on the Ginga light curve. 

To avoid the influence of the accretion flows, we analyze 
only the descending part of the eclipse, which can be con- 
sidered as corresponding to the pure geometrical eclipse 
by the normal component body. 



6. SPECTRAL ANALYSIS 



For spectral analysis, we separated all available observa- 
tions into two parts: below the mean X-ray flux (around 
10 mCrab in the 20-50 keV energy range) and above it. 
The spectra where obtained using the software developed 
by full (IKI RAS). No significant difference in the spec- 
tral shape was found (see Fig. [5)- This suggests a fairly 
broad region emitting in hard X-rays with size compared 
to that of the accretion disk, since (apart from orbital 
eclipses) the observed flux varies due to precession of 
the disk. The spectral shape of this region in the 20- 
100 keV range can be fitted by thermal comptonization 



model (CompPS, fl2lo . Equally good fits were found for 
two cloud geometries: (1) a hemisphere with temperature 
kT = 25 ± 8 keV and optical thickness for electron scat- 
tering T y = 1.8 ± 0.6, and (2) a slab with kT = 23 ± 6 
keV and t v = 2.8 ±0.6 (see Fig. gjand TableGJ. In both 
cases the inclination angle is assumed to be 60 degrees. 
A blackbody component with kT^ = 1.8 keV is added 
to slightly improve the fits. 



This analysis suggests that in addition to thermal X-ray 
emission from jets (which dominates at smaller ener- 
gies in the standard X-ray band, see fl3ll and references 
therein for earlier works), a broad hot corona emitting 
in hard X-ray exists in SS433 around jets, first discov- 
ered by INTEGRAL (Cherepashchuk et al |||, H). From 
the joint analysis of X-ray orbital eclipses and preces- 
sional variability presented in this paper it follows that 
the size of the corona must be comparable with that of 
the accretion disk (10 11 — 10 12 cm). For the physical pa- 
rameters obtained from spectral fits the electron denisty 
around 10 12 cm' 3 can be derived. We note that such a 
denisty is indeed expected in the wind outflowing from a 
supercritical accretion disk with M ~ 10 -4 M Q /yr and 
v ~ 3000 km/s at distances ~ 10 12 cm from the center 
where the photosphere should be formed. The heating of 
such a corona up to 25 keV can be due to violent interac- 
tion of mildly relativistic jets with the surrounding dense 
wind (0). 




Figure 6. The IBIS/ISGRI 10 - 200 keV spectrum of 
SS433 and its best fit by thermal comptonization model 
CompPS with slab geometry (solid line). 
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Figure 7. Eclipse fit for the mass ratio q = 0.3. 



7. ANALYSIS OF ORBITAL AND PRECES- 
SIONAL VARIABILITY 



We carried out the joint analysis of orbital eclipses 
(ingress only) at the precessional phase ip prec ~ 0.1 and 
of the precessional light curve itself. 

For the analysis of X-ray variability we used a geomet- 
rical model of SS433 applied earlier to the analises of 
Ginga data O and the INTEGRAL X-ray light curve 
||5|. A close binary system consists of an opaque "nor- 
mal" star shaped by a Roche equipotential surface and a 
relativistic object surrounded by an optically and geomet- 
rically thick "accretion disk". The size of the "normal" 
star is determined by the sizes of the critical Roche lobe 
for the mass ratio q = M x /M v (here M x is the mass of 
the relativistic object). 



The "accretion disk" includes the disk itself and an ex- 
tended photosphere formed by the outflowing wind. The 
disk is inclined with respect to the orbital plane by the 
angle 20.30°. The opaque disk body with a central cone- 
like region (funnel) is described by the radius r*d and the 
cone opening angle uj. The accretion disk radius is lim- 
ited by the distance from its center to the inner La- 
grangian point. The relativistic object is surrounded by 
a transparent homogeneously emitting spheroid with the 
visible radius rj and height bj, which could be interpreted 
as a "corona" or a "thick jet" (without any relativistic mo- 
tion). 

Only the hot "corona" is assumed to emit hard X-ray flux, 
while the star and the opaque disk eclipse it during orbital 
and precessional motion. Physically, such a "corona" 
could be thought of as hot low dense plasma filling the 
funnel around the relativistic jets. 

During precessional motion the inclination of the disk 
with respect to the observer changes, causing different 
conditions of the "corona" visibility. Observed preces- 
sional variability can thus be used to obtain a "vertical" 
scan of the emitting structure, restricting the height bj 
and the cone opening angle uj. The orbital variability 
(eclipses) allow one to scan the emitting structure "hori- 
zontally", rectricting the possible values of r<2, Vj, uj and 
the mass ratio q. 

The results of the joint analysis of the orbital and preces- 
sional variabilities are the following. 

From the analysis we conclude that in the framework of 
our geometrical model of the hard X-ray emitting region 
in SS433 the eclipse shape (ingress) and the precessional 
amplitude are best reproduced by a wide oblate "corona" 
above an optically thick accretion disk. 

At any q, the best fits for the observed wide orbital eclipse 
d'prec ~ 0.1) are obtained for large visible radii of the 
hot corona rj just slightly smaller than the disk radius 

(rj/r d = 
about 0.3 
70 - 80°. 



0.7 — 0.9). The height of the "corona" bj is 
— 0.6 of rd and the cone opening angle u 



The precessional light curve taken alone does not restrict 
the binary mass ratio q: equally good fits are found for q 
varying from 0.05 to 1.0. 

Analysis of the orbital X-ray eclipse ingress only (ob- 
served at if) P rec ~ 0.1) shows that in the range of q = 
0.1 — 0.6 the fit accuracy varies by some 10% only. 

Simultaneous analysis of both precessional and orbital 
variability of SS433 allows one to constrain q to the range 
0.3 - 0.5. The fit is illustrated in Fig.|7]and Fig.|8] 



8. CONCLUSIONS 



Most probably, SScontains a blach hole. The detection of 
a hot extended corona around the supercritical accretion 





Figure 9. The graphical scheme of the emitting regions 



disk is a new result. 

The future plans are connected with a high resolution 
spectroscopy of SSat the face-on disk phase (planned 
Subaru observations). The new eclipse observations at 
■ijj = — 0.1 in hard X-ray range would confirm the 
growth of the eclipse width with the photon energy, which 
would directly evidence the presence of an extended hot 
corona around the disk. Also the peculiarities of the as- 
cending side of the eclipse are to be studied. These tasks 
are planned to be done by the INTEGRAL observations in 
May 2007. 
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